A MULTI-DIMENSIONAL LIEB-SCHULTZ-MATTIS THEOREM 



BRUNO NACHTERGAELE AND ROBERT SIMS 



o 



Abstract. For a large class of finite-range quantum spin models with half-integer spins, we 
prove that uniqueness of the ground state implies the existence of a low-lying excited state. 
For systems of linear size L, with arbitrary finite dimension, we obtain an upper bound on 
the excitation energy (i.e., the gap above the ground state) of the form (ClogL)/i. This 
«) . result can be regarded as a multi-dimensional Lieb-Schultz-Mattis theorem |14| and provides 

' a rigorous proof of the main result in j8j . 
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1.1. Introduction. Ground state properties of Heisenberg-type antiferromagnets on a va- 
riety of lattices are of great interest in condensed matter physics and material science. An- 
\Q • tiferromagnetic Heisenberg models are directly relevant for the low-temperature behavior of 
^ ■ many materials, most notably the cuprates that exhibit high-T^ superconductivity [16]. 

\ There are several general types of ground states that are known, or expected, to occur 

O \ in specific models: a disordered ground state or spin liquid, critical correlations (power 
^ ■ law decay), dimerization (spin-Peierls states), columnar phases, incommensurate phases, 
. and Neel order. More exotic phenomena such as chiral symmetry breaking have also been 
considered [211 E2]. 

^ \ Which behavior occurs in a given model depends on the lattice, in particular the dimension 

+^ ■ and whether or not the lattice is bipartite, on the type of spin (integer versus half-integer) 
and, of course, also on the interactions. In this paper we are considering a class of half- 
integral spin models (or models where the magnitude of at least some of the spins is half- 
integral). Our aim is to prove a generalization of the Lieb-Schultz-Mattis Theorem |14j . 
^ ' Such a generalization was presented by Hastings in [8] and a substantial part of our proof is 
based on his work. Our main contribution is to provide what we hope is a more transparent 
argument which in addition is mathematically rigorous. 

The well known theorem by Lieb and Mattis [13j implies, among other things, that the 
ground state of the Heisenberg antiferromagnet on a bipartite lattice with isomorphic sub- 
lattices, is non-degenerate. For one-dimensional and quasi-one-dimensional systems of even 
length and with half- integral spin Affleck and Lieb [1], generalizing the original result by 
Lieb, Schultz, and Mattis [H], proved that the gap in the spectrum above the ground state 
is bounded above by constant/ L. A vanishing gap can be expected to lead to a gapless con- 
tinuous spectrum above the ground state in the thermodynamic limit. Such an excitation 
spectrum is generically associated with power-law (as opposed to exponential) decay of cor- 
relations. Aizenman and Nachtergaele proved for the spin-1/2 antiferromagnetic chain that 
if translation invariance is not broken (in particular, when the ground state is unqiue), the 
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spin-spin correlation function can decay no faster than [3J. In other words, uniqueness 
of the ground state imphes slow (power-law) decay of correlations. Recently, it was proved 
rigorously that a non- vanishing spectral gap implies exponential decay of correlations [20119]. 
Therefore, non-exponential decay of correlations implies the absence of a gap. In particular, 
the result by Aizenman and Nachtergaele implies the absence of a gap in the infinite spin- 1/2 
antiferromagnetic chain if the translation invariance is not broken, e.g., if the ground state is 
unique. This result can be generalized to an interesting class of antiferromagnetic chains of 
half-integer spins |T7j. The Lieb-Schultz-Mattis Theorem has also been extended to fermion 
systems on the lattice [251 121] ■ All these results are for one-dimensional systems. The bulk 
of the applications of the spin-1/2 Heisenberg antiferromagnet is in two-dimensional physics 
and therefore, the rigorous proof we provide here, based in part on ideas of Hastings [8], 
should be of considerable interest as it is applicable to higher-dimensional models. 

The most common argument employed to bound a spectral gap from above uses the 
variational principle. Often, the variational state is a perturbation of the ground state. The 
proofs in [2] and [T] are of this kind. However, since the ground state is not known, and no 
assumptions are made about it except for its uniqueness, these proofs are not a variational 
calculation in the usual sense. The variational states are defined by acting with suitable 
local operators A on the (unknown) ground state. 

For a finite volume Hamiltonian Hi generated by a potential $ of the type we consider 
(see the paragraph including (11. 7p and (11.81) in Section [112] for the relevant definitions), and 
with a unique ground state, it is straightforward to show that the gap above the ground 
state, ■jl, is bounded uniformly in L. To see this, note that for any ground state vector Q 
and for any site x, there exists a unitary on the state space of x with vanishing expectation 
in the state Q, i.e., UQ _L Q. Since Q is the unique ground state by assumption, UQ is a 
variational state for the gap. Therefore, we have the bound 

(1.1) < {^,[HL,U]n) <2MY,\MX)\\ <2|||$|||i, 

which is uniform in the system size L. Here, |||$|||i is as defined in (II. 13p . See Section 15.51 
for the proof that such a unitary exists. 

In order to obtain a better bound on the energy of the first excited state one has to exploit 
the few properties assumed of the ground state, such as its uniqueness and symmetries. 
Furthermore, one must show that any proposed variational state has a sufficiently large 
component in the orthogonal complement of the ground state. In Section 12. 2[ we propose a 
variational state for finite systems of size L and then demonstrate the relevant estimates, as 
mentioned above, in Sections [3] and HI It is interesting to note that the energy estimate we 
obtain will itself contain the spectral gap of the finite system in such a way that assuming a 
large gap leads to an upper bound less than the assumed gap. From this contradiction one 
can conclude an upper bound on the finite-volume gaps. 

Our results apply to a rather general class of models, which we will define precisely in the 
next section. The application of our general result to spin-1/2 Hamiltonians with translation 
invariant (or periodic) isotropic finite-range spin-spin interactions on a d-dimensional lattice 
is easy to state. First, let = [1, L] x Vl with L even and with periodic boundary conditions 
in the 1-direction, i.e., in the direction that is of even size. It will be important that the 
number of spins in Vl, \Vl\, is odd, and satisfies iVj,] < cL"^'^, for some d> 1 and a suitable 
constant c. Assuming that the model defined on A^, has a unique ground state, we prove 
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that the spectral gap 7^ satisfies the bound 

,logL 



;i.2) 7L < C- 



L 



where C depends on d and the specifics of the interaction, but not on L. 

Because of the presence of the factor logL, the bound (11.21) apphed to one- dimensional 
models does not fully recover the original Lieb-Schultz-Mattis Theorem in [13] or the bound 
proved by Affleck and Lieb in [1]. This indicates that in general our bound is not optimal. 
Our proof uses in an essential way Lieb- Robinson bounds [IS EDI E], as does Hastings' 
argument in [8] , and the appearance of the factor log L seems to be an inevitable consequence 
of this. In fact, it is known that the standard Heisenberg antiferromagnets with spin > 1 on 
the two-dimensional square lattice or with spin > 1/2 on Z'^ with d > 3, have Neel ordered 
ground states [5], [10] and in that case one can show that the gap is bounded hj C/L (see, 

e.g., mm)- 



1.2. Setup and main result. The arguments we develop below can be applied to a rather 
general class of quantum spin Hamiltonians defined on a large variety of lattices. We believe 
it is useful to present them in a suitably general framework which applies to many interesting 
models. Attempting to be as general as possible, however, would lead us into a morass of 
impenetrable notation. Therefore, we have limited the discussion of further generalizations 
to some brief comments in Section 11.51 

We assume that the Hamiltonians describe interactions between spins that are situated 
at the points of some underlying set A. For simplicity, one may think of A = Z'^, but we 
need only assume that the set A has one direction of translational invariance, which we will 
refer to as the 1-direction. We assume that there is an increasing sequence of sets {A^}^^^ 
which exhaust A of the form A^, = [1,L] xVl where |Vi| < cL'^~^ for some d > 1. Here each 
X G A^ can be written as x = {n, v) where n G {1, 2, ■ ■ ■ , L} and v G Vl, and we will denote 
by (n, Vl) the set of all a; G A^, of the form x = {n, v) for some v G V^. 

Estimates on the decay of correlations in the ground state and Lieb-Robinson bounds on 
the dynamics will play an important role in the proof of the main result. Both are expressed 
in terms of a distance function on A, which we will denote by d. Often, A has the structure 
of a connected graph and d{x, y) is the minimum number of edges in a path from x to y. 
In any case, we will assume that d is a metric and furthermore that there is a function 
F : [0, 00) — >• (0, 00) satisfying the following two conditions. 

Condition Fl: F is uniformly integrable over A in the sense that 
(1.3) := sup^F((i(x,?/)) < 00. 



Condition F2: F satisfies 

ilA) CiF) := sup ^'"'p^''/'^''" < ~- 

x',j/eAf^ F{d{x,y)) 

which means that the "convolution" of F with itself is bounded by a multiple of itself. 

Fl and F2 are restricitve conditions only when A is infinite, however, for finite A, the 
constants and C{F) will be useful in our estimates. It is also important to note that for 



4 



BRUNO NACHTERGAELE AND ROBERT SIMS 



any given set A and function F that satisfies FX and F2 above, we can define a one-parameter 
family of functions, Fx, A > 0, by 

(1.5) F^{x) ■.= e-^^F{x), 

and easily verify that Fl and F2 hold for Fa, with ||Fa|| < ||F|| and Ca(F) < C(F). 

As a concrete example, take A = Z'^ and d{x, y) = \x — y\. In this case, one may take the 
function F{x) = (1 + x)~'^~'^ for any e > 0. Clearly, fll.Sp is satisfied, and a short calculation 
demonstrates that (II ■4p holds with 

Each X G A is assigned a finite-dimensional Hilbert space Tix- For any finite subset X C A, 
the Hilbert space associated with X is the tensor product Hx = ^x^x^^i and the set of 
corresponding observables supported in X is denoted by Ax = B{T-Cx), the bounded linear 
operators over Tix- These local observables form an algebra, and with the natural embedding 
of Axi in Ax2 for any Xi C X2, one can define the C*-algebra of all observables. A, as the 
norm completion of the union of all local observable algebras Ax for finite X C A. Since we 
have assumed that A^^ is of the form [/, r] x Vl with r — / = L — 1, we can define translation 
automorphisms r„, for n G Z, which map A(m,VL) i^^*^ -^{n+myL) for all m G Z. 

An interaction for the system is a map $ from the finite subsets of A to ^ such that for 
each finite X C A, <I>(X)* = G Ax- For given A and F, and any A > 0, let i3A(A) be 

the set of interactions that satisfy 

(1.7) ||<I>||a := sup y J'^^^^\, < 00. 

-•2'eA^^^i^A(rf(x,y)) 

All interactions considered in this paper are assumed to belong to ;Ba(A) for some choice 
of F and A > 0. The constant ||$||a will show up in many estimates. The finite volume 
Hamiltonians are defined in terms of the interaction $ in the usual way by 

(1.8) Hl= $(X) + boundary terms . 

XcAl 

We will always assume periodic boundary conditions in the 1-direction and arbitrary bound- 
ary conditions in the other directions (i.e., any boundary terms in the other directions are 
included in the definition of $). 

The condition that ||$||a is finite is sufficient to guarantee the existence of the dynamics 
in the thermodynamic limit as a one-parameter group of automorphisms on A. In particular 
this means that the limits 

(1.9) af{A) := lim af'^{A) := lim e^'^^Me"**^^ 

L— >oo L— ►oo 

exist in norm for all t G M, and all observables A G Ax, for any finite X C A. We will often 
suppress the L or $ dependence in the notation af See [H [231 HB] for more details. 

Next, we turn to a set of conditions that more specifically describe the class of models to 
which the Lieb-Schultz-Mattis Theorem may be applied. 

Condition LSMl: We assume that the interaction is translation invariant in at least one 
direction, which we will take to be the 1-direction. This means 

(1.10) <l>(X + ei)=ri(<l>(X)), 
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where, for any X C A, X + ci is translation of all points in X by one unit in the 1-direction. 
We will consider finite systems with Hamiltonians defined with periodic boundary con- 
ditions in the 1-direction. For convenience of the presentation we will assume free boundary 
conditions in the other directions but this is not crucial. Since we have assumed periodicity 
in the 1-direction, we can implement the translation invariance for finite systems by a unitary 
T e Aa^ such that $(X + d) = T*$(X)r, for all X C A^. Here T depends on L, but we 
suppress this dependence in the notation. 

Condition LSM2: The interactions are assumed to be of finite range in the 1-direction, i.e., 
there exists R> (the range), such that if X C A and X 3 Xi = {ni,Vi) for i = 1,2 with 
\ni - nsl >R, then $(X) = 0. 

Condition LSM3: We assume rotation invariance about one axis. More precisely, we assume 
that there is a hermitian matrix in every A^x}, a; G A, which we will denote by S^, with 
eigenvalues that are either all integer or all half- integer (i.e. belonging to Z + 1/2). We also 
require that Tm{S^) = S^+mei- Define, for G M, the unitary U{9) G Aaj^ by 

(1.11) U{e) = (g) e''^' . 

The interaction is taken to be rotation invariant in the sense that for each finite X C A 

(1.12) U*{e)^{X)U{e) = <I>(X) for all G M . 

Condition LSM^: We assume that the are uniformly bounded: there exists S such that 
< for all x G A. The following condition, which we call odd parity, is crucial: define 
the parity of x, px to be if the eigenvalues of S^. are integers, and px = 1/2 if they are 
half-integers. We assume that ^^gvi Pin,v) G Z -|- 1/2, for all n G Z. The simplest and most 
important case where this is satisfied is when we have a spin 1/2 at each site, and \Vl\ is 
odd. 

Condition LSM5: The ground state of Hi is assumed non-degenerate. This implies it is 
an eigenvector of the translation T and rotations U{6). Without loss of generality we can 
assume that 1 is the corresponding eigenvalue of T (if the eigenvalue is e*"'^, replace T by 
e-i<^T). We also assume that the ground state has eigenvalue 1 for the rotations U{6). 

Condition LSM6: We assume that there are orthonormal bases of the Hilbert spaces 
with respect to which 5*^ and $(X) are real, for all x G Al,X C A^,. This condition is only 
used in the proof of Lemma [2. 1[ Therefore, this condition may be replaced by the property 
proved in that lemma. 

We will also use the following quantities: 

(1.13) |||$|||i := sup^||<l>(X)|| < oo, 
and 

(1.14) |||$|||2 := supJ]|X| J]1|[5,3„<|.(X)]|| <oo. 

It is not hard to show that the conditions Fl and F2 are sufficient to imply that |||$|||i and 
1 1 1$| 1 12 are finite. 
We can now state our main result. 
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Theorem 1.1. Let 7/, be the spectral gap, i.e., the difference between the lowest and next- 
lowest eigenvalue of the Hamiltonian of a model satisfying conditions Fl, F2, and LSMl- 
6. Then, there exists a constant C , depending only on properties of A ( such as the dimen- 
sion), the constants \\F\\ and C{F), and the interaction (\\^\\\, for some A > 0, HI^IHi, and 
1 1 1"^! I such that 



1.3. Structure of the proof. The simplest way to present the proof IS cLS cL proof by 
contradiction. Under the assumption that there exists a sufficiently large constant C > 0, 
such that 'Jl exceeds {C log L)/L for large L, we will construct a state orthogonal to the 
ground state with an energy difference that is boundable by a quantity that is strictly less 
than the assumed gap for sufficiently large L. Thus, the proof is in essence a variational 
argument. The variational state is constructed as a perturbation of the ground state, as the 
solution of the differential equation proposed by Hastings [8J with the ground state as initial 
condition (see Section [2l in particular fl2.25p . for this equation). The important idea is that 
this equation will lead to a state which resembles the ground state of the Hamiltonian with 
twisted rather than periodic boundary conditions (see Section 12.11 for the definition of the 
twists), at least in part of the system, say the left half. In the right half the ground state 
will be left essentially unperturbed. This state is defined in Section [2l 

After the variational state has been defined, there are two main steps in the proof: estimat- 
ing its excitation energy and verifying that it is "sufficiently orthogonal" to the ground state. 
In general, one may also have to consider the normalization of the variational state, but in 
our case the differential equation defining it will be manifestly norm preserving. Hence, this 
is not an issue for our proof. 

The main difficulty is that under the general assumptions we have made, no explicit 
information about the ground state is available. Its uniqueness, translation, and rotation 
invariance are the only properties we can use. In combination with the general assumptions 
on the interactions and the assumption on the magnitude of the spectral gap above the 
ground state, however, one can obtain an upper bound on the decay of correlations of the 
ground state in the 1-direction. The recently proved Lieb-Robinson bounds [201 El HB] will 
be essential to show that the effects of the perturbations we define in the left half of the 
system remain essentially localized there. This allows us to compare the energy of the 
variational state with the ground state energy of a Hamiltonian, introduced in (12. 5p . 

which, instead of twisted boundary conditions, has two twists that cancel each other. The 
twisted Hamiltonian is unitarily equivalent to the original one and therefore has the same 
ground state energy. We work out this argument in Section [31 The result is 



where ijji is the normalized variational state we construct, and Eq is the ground state energy. 
The dependence of both quantities on L is surpressed in the notation, i/, C and c are positive 
constants that only depend on properties of the lattice and the interactions. The correction 
terms appearing above, and also in (I1.17P below, can be made explicit by the estimates 
provided in Section [51 They depend on the quantity 7lL in such a way that assuming there 
exists a constant C > for which 'JlL > C for sufficiently large L, they are uniformly 
bounded in L. Due to the nature of our proof of Theorem II. H see below, we do not write 
these additional terms out explicitly. 




1L<C 



logL 



L 



(1.16) 
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For the orthogonality, our strategy is to show that ipi is almost an eigenvector of the 
translation T with eigenvalue —1. Since the ground state ipo is an eigenvector of T with 
eigenvalue 1, by assumption, this shows that ipi is nearly orthogonal to ipo- In Section H] we 
obtain a bound on their inner product of the form: 

(1.17) |(V^i,V^o)| < C'L^'e-^'""^^ (1 + corrections) , 

where u', C and c' are positive constants similar to C and c. 
The proof of Theorem 11.11 then easily follows. 

Proof, (of Theorem 11.11 ) Suppose that '-y^L > ClogL with a sufficiently large constant C. 
In this case, the correction terms which appear in the bounds (11.161) and (I1.17P above are 
negligible. It is easy to see then that one obtains a contradiction for L large enough. □ 

To help the reader see the forest through the trees we have tried to streamline the estimates 
in Sections [3] and H] by collecting some results of a more technical nature in Section [51 

1.4. Examples. The conditions LSMl-6 we have imposed on the models are not unrea- 
sonable. We will illustrate this by considering various antiferromagnetic Heisenberg models 
defined on Al = [1, L] x Vl, where for each L, Vl is a. finite set. As before, at each x G A^, 
we have a finite-dimensional spin with spin- matrices 5* , i = 1, 2, 3, and we consider Hamil- 
tonians of the form 

(1.18) H= Yl J{x,y)S^-Sy, 

where J(x, y) G M are the coupling constants. 

If Vl C Z'^~^, with d > 1, and such \Vl\ < cL'^~^, for a suitable constant c, which 
describes the case for (i-dimensional systems defined on subsets of Z*^, there exists a function 
F satisfying Conditions Fl and F2 as we have indicated in the paragraph containing (II. 6p . 
It is also easy to see that if Vl is a fixed finite set independent of L, in which case the system 
is (quasi) one-dimensional, any function F that works for the one- dimensional lattice will 
suffice. All the examples we discuss below will be of this type. 

Certainly, there are still many Hamiltonians of the form (11.181) that fail to satisfy all 
six conditions, but this is generally for a good reason. For example, without translation 
invariance in at least one direction one can easily have a non- vanishing gap above the ground 
state. 

Condition LSM2, finite-range, does not need to be satisfied in the strict sense. Sufficiently 
rapidly decaying interactions could also be treated. For the present discussion, let's assume 
that the model is translation invariant in the 1-direction, and that the interactions are 
nearest neighbor in the 1-direction in the sense that for any x = [rii, v) and y = (n2, u), with 
\ni — > 1, we have J(x, y) = 0. 

The rotation invariance about at least one axis imposed in Condition LSM3 is essential for 
the type of result we prove. The models (11.181) have full rotation invariance, so they clearly 
satisfy this condition. Anisotropic models of the XXZ type would still satisfy LSM3. 

In order to satisfy LSM4, we have to assume a uniform bound on the size of the spin. 
Clearly, all models with only one kind of spins or a periodic arrangement of spin magnitudes 
satisfy this condition. Since we already assumed translation invariance in the 1-direction, 
we can verify the odd parity condition by adding the magnitudes of all spins in the "slice" 
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(1,Vl). If we have only spin 1/2's, e.g., we simply need that \Vl\ is odd. For the one- 
dimensional chains of identical spins of magnitude S, the condition requires that S is half- 
integral. Haldane's Conjecture [U [7] predicts that for integer values of S there exists a 
non- vanishing gap. There are examples of isotropic integer-spin chains which satisfy all the 
other conditions and for which the existence of a non-vanishing gap has been rigorously 
established, such as the AKLT chain [2]. For p-periodic spin chains with a repeating pattern 
of spin magnitudes Si, . . . , Sp, LSM4 is satisfied if S'l + ■■■ + S'p is half-integral. Similarly, 
for spin ladders LSM4 is satisfied if the total spin in each rung is half-integral. 

There is a large class of models for which the uniqueness of the ground state demanded 
by LSM5 follows from the Lieb-Mattis Theorem [13]. For Hamiltonians of the form fll.lSp . 
a simple case where the Lieb-Mattis Theorem applies is the following: if A/, is the union of 
two disjoint subsets A^ yi and A^ ^ of equal size, with J{x,y) < whenever x and y do not 
belong to the same subset, and sufficiently many J{x,y) are non-vanishing such that the 
graph formed by the edges with non-zero coupling constants is connected. This is satisfied if 
Vl C 7a'^~^ is connected and the Hamiltonian is the usual nearest neighbor antiferromagnetic 
Heisenberg model. 

All models of the form (fTTSj) satisfy LSM6. 

The above discussion demonstrates that there is a large variety of models that satisfy 
all conditions of our main theorem. In particular, all nearest-neighbor half-integer spin 
Heisenberg antiferromagnets defined on subsets A^ = [1, 2L] x V;, of d- dimensional hypercubic 
lattice with \Vl\ odd and such that \Vl\ < cL^, for some a > (it is natural but not 
necessary to assume a = d—1), have a unique ground state with a gap 7^ above it satisfying 
7i < C(logL)/L, for some constant C. 

1.5. Generalizations. One can envision several generalizations of Theorem 11.11 An obvious 
one is to remove the condition that the interaction is strictly finite range in the 1-direction. 
It is not hard to see that the arguments given in the following sections can be extended to 
long-range interactions with sufficiently fast decay. 

One may wonder whether the assumption that L is even is essential. It is used in the proof 
of near orthogonality of the variational state, which is based on investigating the behavior 
under translations of the state: the variational state is close to an eigenvector with eigenvalue 
— 1 of the translation operator T, whereas the ground state has eigenvalue 1. Our proof of 
this fact assumes that the ground state is an eigenvector of the rotations with eigenvalue 
1. For L odd, our assumptions preclude the existence of such an eigenvector. However, it 
seems plausible that for odd L a slight modification of our proof will work to show that the 
ground state and the variational excited state have opposite eigenvalues for translations. 

The main applications we think of are to SU(2)-invariant Hamiltonians with antiferromag- 
netic interactions. Affleck and Lieb |lj pointed out that their proof easily extends to a class 
of models with SU(N) symmetry. There are no obstructions to generalizing our arguments 
to such models with SU(N) symmetry given by suitable representations. 

It may also be of interest to consider different topologies of the underlying lattice and/or 
the twistings. Instead of cylindrical systems with periodic boundary conditions which can 
be deformed by a twist, one could apply a similar strategy to systems defined on a ball or a 
sphere. We do not explore such possibilities here. 

Another question we do not address in this paper is under what circumstances the trial 
state we construct is actually a good approximation of a low-lying eigenstate with energy 
close to the first excited state, or even whether it is a state orthogonal to the ground state 
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and with energy bounded by C (log L) / L. We do not believe that statements of this kind hold 
under the general conditions we impose. It is expected that in some cases the true gap of 
the system is much smaller than the bound we prove. This is of course not in contradiction 
with our result, but under such circumstances our method does not provably construct a 
good variational state. There is no reason to assume that it always would. 

2. Construction of the variational state 

2.1. Twisted Hamiltonians. The main motivation behind the construction of the varia- 
tional excited state is that it should resemble the ground state of the model with twisted 
(as opposed to periodic) boundary conditions. Therefore, we first describe some elemen- 
tary properties of a family of perturbations of the Hamiltonian, which we will call twisted 
Hamiltonians for reasons that will become obvious. 

Given an interaction $ which satisfies the general assumptions outlined in Section II. 2[ 
we will now define a two parameter family of "twisted Hamiltonians" to analyze. These 
Hamiltonians will be defined on a finite volume A/, = x V^^ where [1,/^] is considered 
with periodic boundary conditions for some even L > AR, where i? > is the range of 
$ in the 1-direction. Let $l be the periodic extension of $ restricted to A/,. Recall that 
each point x E Al can be written as x = {n, v) where n G {1, 2, . . . , L} and v G Vl, and 
we will denote by (n, V/^) = {a; G Al : x = {n,v) for some v G Vl}. For any ^ G M and 
n G {1,2,..., L}, define the "column" rotations Un{0) by 

(2.1) U4e)= (g) e^'^'. 

xe{n,VL) 

For mG {1,2,. ..,L — 1}, we will denote by Vrn{6) the unitary given by 

(2.2) Vm{e) = (g) Un{e). 

m<n<L 

The "twisted Hamiltonians" are defined to be perturbations of the initial Hamiltonian 
with periodic boundary conditions defined by 

(2.3) H = J2 ^LiX). 

XcAl 

The perturbations have the following form: 

(2.4) He{m) := V;„(^^)*<l>(X)V;,(e) - $(X). 

XcAl 

for m G [R,L — R] to avoid interactions across the seam created by identifying L + 1 with 1. 
Note that here we use the original potential $, and not its periodic extension Clearly, 
if ^ C U^<„<L(n, Vl) or X C Ui<„<^(^, Vl), then VUO)*HX)VUO) - $ (X) will vanish 
by rotation invariance of the interaction, and therefore only those interactions across the 
column (m, Vl) contribute in fl2^ . For 6,6' eR and me {R,R + 1, . . . , L/2 - R} fixed, we 
define 

(2.5) He,e' := H + He{m) + He'{m + L/2), 

to be a doubly twisted Hamiltonian. With m fixed, we regard A^ as the disjoint union of 
two sets 

(2.6) AL = AruAf, 
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where A^^"* consists of two windows, one about each column at which a twist has been 
apphed; namely 

(2.7) Kf^ := Kf\m)UKf\m + L/2) ^n<\ Kf\y) := (J 

\n'y\<T-R 

for y G {m,m + L/2}. Moreover, comprises the remaining strips in Al. Given this 
decomposition of the underlying space, the twisted Hamiltonian can be written as 

(2.8) He,e' = Hj,^} + H^'\ 
where 

(2.9) H^'^= Yl "^LiX), 

and Hq^} denotes the remaining terms in He gi which, due to (12.91) . are supported strictly 
within the windows. 

There are a variety of useful symmetries the Hamiltonians He^gr, introduced in ( 12. Sp . 
possess. With m G {R, R + 1, . . . , L/2 — R} fixed as above, one may define 

(2.10) W{<P) := (g) [/„(-0), 

m<n<m+L/2 

for any real (p. See (12. ip for the definition of the column rotations It is easy to check 
that for any angles 6, 6', G M, one has that 

(2.11) W*{<j))He,e'W{<j)) = He-^,e'+^, 

due to the (term by term) rotation invariance of the interactions. Given this relation, it 
is clear that along the path 9' = —9 the twisted Hamiltonian is unitarily equivalent to the 
untwisted Hamiltonian, i.e., 

(2.12) W{9yHe,-eW{9) = Ho,o = H, 

which, due to the periodic boundary conditions, is not true for general pairs 9, 9'. 

The untwisted Hamiltonian is translation invariant (in the 1-direction), and it is important 
that the twisted Hamiltonians inherit a "twisted" translation invariance. Define 

(2.13) Te,e' = TUU0)Um+L/2{9'). 

where T is the unitary implementing the translation by 1 in the 1-direction. It is then 
straighforward to verify that 

(2.14) Hg^0i = Tq q, Hq^iTq^qi. 

Note that under the odd parity condition LSM4 we have 

(2.15) T2,,o = -T, 

which will be important in the proof of the almost orthogonality of the trial state in Section 
il 

If we denote by ipo the (unique) ground state of if, i.e., Hipo = Eq^Pq, then by translation 
invariance, and specifically LSM5, we have that TipQ = ipQ. Moreover, using the unitary 
equivalence (12.121) . we see that the ground state of the twisted Hamiltonian Hq _0 satis- 
fies He^-eM^.-e) = Eo{9, -9)M0, -0) with Eo{9,-9) = Eo and M0,-0) = V(^)^o. 
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Although the twisted ground state ipoi^^^ ~d) is not translation invariant, it does satisfy in- 
variance with respect to the twisted translations, i.e., Tq_0iPq{9,—9) = iIjq{6,—6). As a 
consequence, we have the following simple but important property of Eq. 

Lemma 2.1. Let Eq{6,6') denote the ground state energy of Hqqi .Then, the partial deriva- 
tives of Eq vanish on the line 9' = —6: 



Proof. First, we note that Eq is differentiable in its two variables in a neighborhood of (0, 0) 
by the non-degeneracy condition LSM5. By unitary equivalence Eq is then differentiable in 
a neighborhood of the line {9, —6). For G M, let £{'ip,(j)) = Eq{iP — cl),ip + (p) denote 
the ground state energy of H^_(j,^^+^. Due to the unitary equivalence eq. fl2.1ip . £ depends 
only on ip. Hence, d^S{ip,(j)) = 0, for all ip,4>. Under the additional assumption that the 
interactions are real (LSM6), we have that Hgg/ = H^g^^g', and therefore £{ip,0) = 

S{—ip,0). Hence, S is an even function of and d^S^ip, (p)\ip=o = 0. Using these properties 
and the fact that the partial derivatives of Eq are linear combinations of the partial derivatives 
of S, we find that both partial derivatives of Eq vanish on the line 6' = —9. □ 

2.2. The variational state. Our aim is to construct a state that resembles the ground 
state of Hg _g in a region surrounding those spins that were twisted by an angle of 6, while 
it otherwise resembles the ground state of H = Hq q. 

From the unitary equivalence fl2.12p we have that Eq{6, —6) is independent of 6*, i.e., 
dgEQ{6, —9) = 0. Moreover, the partial derivatives of Eq vanish on the line {6, —0), as was 
proven in Lemma 12.11 This property, in general, allows one to derive an equation for the 
ground state. 

Consider a differentiable one-parameter family of self-adjoint operators H{x), x G [a, b] C 
M, and let Eq{x) denote the ground state energy of H{x) with a differentiable family of 
ground state eigenvectors ipoi^)- Suppose OxEq^x) = for x G [a, 6]. Then, it is easy to see 
that 'ipQ^x) ± {d^H {x))iPq{x) , from which we obtain: 



where PJ is the spectral resolution for H{x) and af^ is the imaginary-time evolution corre- 
sponding to the Hamiltonian. 

Motivated by this calculation, we introduce the family of operators B{A, H), where H is 
a Hamiltonian for which the dynamics {at \ t G M} exists as a strongly continuous group of 
^-automorphisms and A is any local observable, defined by 



(2.16) 



d,EQ{9,-9)=d2EQ{9,-9)=0. 




[Hix)-EQix)]-' d,H{x)^Q{x). 



(2.18) 




(2.19) 
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where is the imaginary time evolution generated by H. For unbounded Hamiltonians H, 
it may not be obvious that B{A, H) can be defined on a dense domain. However, if ■?/' is a 
ground state corresponding to the Hamihonian H, then B{A, H)ip exists. Moreover, from 
fl2.18p . we conclude that 



(2.20) d,M^) = B{x)M^), 

where B{x) = B {dxH{x), H{x)). Similarly, in the density matrix formalism, for 

(2.21) po{x) := \Mx)){Mx)\, 
equation (12.201) implies that 

(2.22) d^poix) = B{x)poix) + po{x)B{xy 



We will define the proposed excited state ip as the solution of a differential equation anal- 
ogous to (12.201) . First, we need to introduce some further notation. Let H he a Hamiltonian 
for which the dynamics {at} exists; finite volume is sufficient. For any a > 0, t G ]R\{0}, 
and local observable A & A, we may define 

(2.23) A,{tt,H) := —e-^' j aM)j—-^ds. 
In addition, for T > the quantity 

(2.24) BaAAH):=-[ A^iit , H) - A^iit , H)* dt , 

Jo 

will play a crucial role. In Lemma [5.121 of Section O we will show that when projected onto 
the ground state of a gapped Hamiltonians H, the quantity Ba^riA, H) well approximates 
B{A,H) for a judicious choice of parameters, e.g., a = 'Jl/L and T = L/2; we note that 
the observable A must also satisfy the constraint that its range is orthogonal to the ground 
state. With this in mind, consider the solution of the differential equation introduced by 
Hastings in [8]: 

(2.25) dgllJa,T{0) = BaAG)^aA0)^ 

where Ba,T{6) = Ba^xiddHg^Q, Hg^^g), subject to the boundary condition ipa,T{0) = V'o- Note 
that Ba,T{0) is anti-hermitian, and hence any ipa,T{d) solving (I2.25P will have constant norm. 

To be explicit, the proposed state ipa^riO) differs from the actual ground state iIjq{6, —9), 
of the doubly twisted Hamiltonian Hg _g, in three essential ways. Compare (12.191) in the case 
that A = dgHg _g and H = Hg^^g with (I2.24p given that A = dgHg Q and H = Hg _g. 

i) We have introduced a cut-off at T < cxd. 

ii) We have approximated the imaginary-time evolution of an observable A, ait{A), by 

Aa{lt,H)-Aa{zt,Hy. 

iii) We have replaced the observable dgHg^^g with dgHgfl. 

The modifications i) and ii) are of a technical nature, i.e., to make the relevant quantities 
well-defined and amenable to estimation (see Section [5]). The motivation behind the third 
replacement is an attempt to approximate the ground state of the singly twisted Hamiltonian 
He,o- 
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3. Energy estimates 

As is discussed in the introduction, the goal of this section is to prove an estimate of the 
form 

(3.1) \{^,,HL^^)-Eo\<CL''e-''rLL^ 

see (11.161) and Theorem 13. 31 below, for the proposed variational state. Explicitly, we will take 
i^i = '?/'a,T (2vr), i.e., the solution of (12.251) evaluated at ^ = 27r, with the specific choice of 
parameters a = 'Jl/L and T = L/2. Since the operator Ba^riO), defined in (12.251) . is anti- 
hermitian, it is clear that ipi remains normalized, and the bound stated above demonstrates 
that if the gap is sufficiently large, 7^ > Clog(L)/L, then ipi corresponds to a state with 
small (depending on C) excitation energy. An estimate of the form (13.11) . with correction 
terms, can be proven based on the general results in Section [51 In the proof of Theorem 
II. H which is a proof by contradiction, we will be assuming 7^, > Clog(L)/L. Therefore, 
we can assume here, without loss of generality, that there exists a constant c > such that 
'JlL > c for sufficiently large L. This assumption, which is not necessary, will simplify the 
presentation in Sections [3] and HI 

3.1. Local estimates on the states. In this subsection we prove a technical result which 
estimates, uniformly in 6, the norm difference between the ground state of Hg _g and the 
proposed state in the left half of the system, more precisely in the window centered around 
the location, (m, Vl) where the ^-twist occurs. Since the restrictions of the states to the 
half-systems are described by density matrices, it is natural to use the trace norm for this 
estimate. Recall that for any bounded operator A on a Hilbert space 7i, the trace norm is 
defined by 

(3.2) Pill =Tr 71^4, 

assuming this quantity is finite. Using the polar decomposition for bounded linear operators, 
it is easy to see that, alternatively, 

(3.3) Plli= sup \TtAB\, 

\\B\\=1 

Recall that the density matrix corresponding to the ground state of the Hg^^e Hamiltonian 
satisfies the equation 

(3.4) dgpo{6, -9) = B{9) po{9, -9) + p,{9, -9) B{9)\ 

compare with (I2.22p . where we have used the notation B{9) = B{A{9), Hg^^g) for the operator 
B{A, H) as defined in ( 12.19^ and the observable A{9) = dgHg^^g. Note that by construction 
Pq{9, —9) remains normalized. We will often write A{9) = ^1(6') — A2{9) where, the observ- 
ables Ai{9) = dgHg^o and ^2(6*) = dgH^^g are supported in the window about the twists of 
angle 9 and —9, respectively. Regarding Hg gi as a function of two variables, we may write 
Ai{9) = diHg ^g for convenience. The notation Bi{9) = B{Ai{9), Hg _g) will also be useful. 
The proposed state is the solution of 

(3.5) dgPaAS) = [BaA0)^PaA0)]^ 

where the operator Ba^T{9) = Ba^T{Ai{9), Hg^^g) as in (12.241) with observable Ai{9) = 
diHg^-g. The parametrization we choose is a = 7l/L and T = L/2. Since the opera- 
tor Ba^T{9) is anti-hermitian, the solution Pa,T{9) is a density matrix. We will denoted by 
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Trm':[-] the partial trace over the Hilbert space corresponding to A^^'' U A^^'*(m + L/2). Note 
that the terms in the Hamiltonian that have been twisted by an angle 9 are supported in the 
complementary region K^^\m). Given a gap 7^ > above the ground state of the H = ifo,o 
Hamiltonian, we will be able to estimate the trace norm of the difference in the two states 
restricted to A^^\m). We will show this by estimating deTim'^ [pa,T{9) — Po{0, —9)]. 

Theorem 3.1. As described in the introduction, we assume Fl, F2, and LSMl-6. If there 
exists a constant c > such that 'JlL > c for sufficiently large L and we choose the parameters 
a = 'Jl/L and T = L/2, then there exists constants C > and k > so that 

(3.6) sup \\TT„ApaA0)-po{d,-O)]\\, < CL^'^e-''^^'^, 

e6[0,27r] 

for L large enough. Here C and k depend only on the interaction $ and the underlying set 
A. 

We note that the assumption concerning the existence of a constant c > such that 'JlL > 
c for sufficiently large L is not necessary. We impose it here for simplicity of presentation. 
Without this additional assumption, one may prove an analogue of (13.61) . which contains 
correction terms, by inserting the bounds proven in Section [5] directly into the proof given 
below. Since we make this assumption, it is convenient to state a lemma which compiles many 
of the technical results found in Section[5]and applies them to the present set-up. For this, we 

need two more definitions. Denote by Bj^\9) the operator defined by Ba^x (^■Ai{9) , Hg^l 

where the Hamiltonian Hg^^ is the full Hamiltonian Hg^^g restricted to the windows about 
the twists, see (12. 8p . Lastly, set Pq to be the projection onto the ground state 4'o{9, —9) of 
the twisted Hamiltonian Hg^^g. 

Lemma 3.2. Assume Fl, F2, and LSMl-6. If there exists a constant c > such that 
IlL > c for sufficiently large L and we choose the parameters a = •yi/L and T = L/2, then 
there exists constants C > and k > for which both 

(3.7) sup \\B,A9) - B;^^\9)\\ < CL^'^e-^^-^ 

6»e[0,27r] 

and 

(3.8) sup UBa,T{9) - B,{9)) P^W < CL'^e"^^^^ 

6'e[0,27r] 

when L is large enough. 

Proof, (of Lemma 13. 2p Equation (13. 7p follows by combining Lemma 15.31 and Remark 15.71 
Using Lemma [5.121 and Remark 15. 13^ one obtains (13.80 . □ 

Proof, (of Theorem 13. ip The proof of Theorem 13.11 follows by deriving a uniform bound on 
the 6'-derivative of the differences in these density matrices. Specifically, the bound is in 
trace norm, and the uniformity is with respect to ^ G [0, 27i]. 

Using (13. 4p . (13. 5p . and inserting the local operator B'^\9) for comparison, one may easily 
verify that 

3 

(3.9) dgTTm4paA^)~Po{0,-0)] = Tr„c Q (0) , p,,^ (^) " Po (^, -^^) ] ) + $^r,(^), 

i=l 
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(3.10) 



r^ie) := Tr„c ( Po(^,-^)] - 9ipo(^,-^)) 



(3.11) 
and 

(3.12) r3(0) := Tr^e [d^p,{e, -9) - dep^{e, -9)] . 

As Ai{6) is supported near (m, V/^) and Hg^l contains only those interaction terms over 
sets X C A^^"*, it is clear that Bj^\6) is contained in the algebra of local observables with 
support in A^^\m); we will denote this algebra by A{m). Therefore, 



(3.13) Tr„c ( B;!P{e),p,A0) - Po{0, -i 



i?iy(e),Tr„c {p^A(^)-p,{e,-e)) 



Since B^\6) is anti-hermitian, we may apply norm preservation, i.e. Theorem 15.151 to 
(13. 9p and conclude that 

(3.14) iiTr^^[p^^^^e)-poie,-e)]\\^ < V f \\r,ie')hde'. 

i=i 

We need only bound the trace norms of the remainder terms ri{6). 

As Pa,T{d) is a density matrix, in particular non-negative with a normalized trace, one has 
that 

(3.15) Wnmi < 2||5.,T(^)-5iy(^)|| < CL2'^e-^^-^ 

using Lemma [3.21 above. 

To estimate r2{0), we note that as in (13. 4p . 

(3.16) d^po{e,-e) = Bi{e)po{e,-e) + po{e,-e)Bi{e)\ 

where di denotes differentiation with respect to only the first twist angle, namely 6, which is 
situated near the sites (m, Vl)- Here we have also used that diEo{6, —6) = 0, see Lemma 12111 
A simple norm estimate yields that 



r2 



(3.17) 



1 < 2 
< 2 



BaAO) - B!S^\e) + 2 \\{B^AO) - BmP^ II • 



Appealing again to Lemma [3.21 we see that r2{6) satisfies the desired bound. 

Lastly, r3 (6*) = Tim': [d2Po{0, —9)]. Since we have shown in Lemma 12. II that doEn(9. —i 
as well, the analogue of (I3.16P holds for d2Po{6, —6). Thus, 



^3 



sup |Tr[o {B2{e)po{e,-e) + po{e,-e)B2{ey)]\ 



l|0|| = l 



(3.18) 



< 2 sup 



( ^o(^^, -e), O aa {A2{9)) ^o(^^, -9) ) \ dt, 



OeA{m): Jo 

IIO||=l 



where the observables O are arbitrary elements of ^(m), again, the algebra of local ob- 
servables with support in A^^\m). Integrals of this type are bounded using Lemma [5^ 
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see also Remark I5.10[ Since the observables we are considering have a separation distance 
proportional to L, we may estimate 

(3.19) Iks Will < CL^'^e-'^^^^. 

Combining the results found on each of the remainders, we arrive at the estimate claimed 
in ([MD. □ 

3.2. Bound on the energy. Equipped with Theorem 13.11 and Lemma 13.21 we may now 
bound the excitation energy corresponding to the proposed state. 

Theorem 3.3. Assume Fl, F2, and LSMl-6. If there exists a constant c > such that 
IlL > c for sufficiently large L and we choose the parameters a = '-jl/ L and T = L/2, then 
there exists constants C > and k > so that 

(3.20) |(^i,i^L^i) - ^ol < CL^'-^e-^^^^ 
for large enough L. Here, we take ipi = il)a,T{'^'^)- 

The proof of this theorem may be understood as follows. Recall that the ground state 
energy of the doubly twisted Hamiltonian is independent of 6, i.e., 

(3.21) Eo = (^o,^L^o) = {Md.-0),He,-eMd,-0)). 

Moreover, the separation between the twists of angle 9 and —9 grows with the volume. 
Locality should enable one to estimate the energy difference between performing two twists, 
the ground state, and performing only one twist, the excited state. A rigorous version of 
this idea is described below. 

First, we recall some of the notation introduced in Section 12. 1[ We have written the 
twisted Hamiltonian as the sum of two terms 

(3.22) He,-e = H^, + if(^). 
It is useful to further subdivide the twisted terms as 

(3.23) Hi^_l = HPim) + H^^\m + L/2), 

where H^^^ (m) contains all those interaction terms in Hg^l with support in a window about 

the twist of angle 9, i.e. A^^\m), and similarly, H^\m + L/2) contains all those interaction 

terms in Hg^l with support in A^^\m + L/2). The untwisted terms in fl3.22p are supported 
in the remaining strips. We refer to equations (12. 5p - (12.91) for more details. It was also noted 
in Section HIT] that 

(3.24) W{9y He,-eW{9) = H. 
see (EH. 

Now, for any state ip, one may calculate the expected energy due to a single twist: 

(3.25) = ^0 + RiiO) + R2{9). 

Here, we inserted appropriate terms so that we may compare {ip, Hg^ip) to the ground state 
energy; the remainder terms are given by 

(3.26) := {^,H^^\m)i,) - {^0, -0), im)MO, -0)) 
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and 



(3.27) R2{e) := (i^rV + L/2) + H^'^) ^) - (^o, {H^^\m + L/2) + H^'^) ^o), 
where iIjq is the ground state of H = Hq q. The bound 



(3.28) 



{ij,He,oi^) - Eo \ < \Ri{9)\ + |i?2(^^)|, 



readily follows for any state ip. 

Proof, (of Theorem 13.31) For each fixed 6, the bound (13.281) is valid for our proposed state 
'^a,T{S)- We will estimate the resulting remainders uniformly for 6 G [0, 27r] and thereby 
prove the claimed result. 

To see this, we first rewrite the remainders in terms of the density matrices of the states 
restricted to the region containing the first twist. It is clear that 



(3.29) 



RiiO) 



Tr 
Tr, 



Ti^4PaA0)-PoiO,-O)] 



m] 



where the partial traces defined just prior to Theorem 13. 1[ Thus 



(3.30) 



\Rim < ||if(^)(m)i|||Tr^c[p,,r(^)-po( 



where we have used Theorem 13.11 In fact, from the assumptions we have made, one verifies 
that 



(3.3i; 



\Hr\m) 



< 2 ^||<I>(X)|| < 2|||<l>|||i|Af'^(m)| < CL'^-\ 



For the second remainder. 



(3.32) 



R2{e) 



Tr 



(Pa,T(^)-po(0,0)) (H, 



L/2) + if(^) 



we observe that the only 6 dependence is in the density matrix corresponding to the proposed 
state. Using the differential equation (13.51) . we find that 

(3.33) R'M = i:r{^[BaA0).Pa,Tm[H^^\rn + L/2) + H^') 



= - Tr 

The first term above is easy to estimate. Recall that the quantity b'^\6) is supported in 
A^^\m), whereas H^\m + L/2) has support in A!"^\m + L/2). Thus 



(3.34) 
and moreover, 
(3.35) 

< 2 



B,A(^),H^^\m + L/2) = BaA0) - BrT'i0),H'J'>{m + L/2) 



Tr 



B,A(^) - B^^\e),H^''\m + L/2) p^^^) 



H^^\m + L/2) B^AO) - B':'\e) 



{W)i 
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The second term may be similarly estimated. Let Hg^l be defined as in (12. 8p . excepting 
that the windows are slightly smaller: of size j — 2R. Then [B^^\e), H^^^] = 0, and the 
argument above applies. We have bounded R2{d)- D 

4. Orthogonality 

We will now prove that, under the assumptions given in the introduction, the proposed 
state is nearly orthogonal to the ground state. As in Section El we again make the assumption 
that > c > 0, for sufficiently large L. 

The reasoning behind orthogonality is simple. From LSM5, we know that the ground state 
is an eigenvector of the translation operator with eigenvalue 1, i.e., TipQ = ipQ. On the other 
hand, the proposed state will very nearly be an eigenvector of T2Ttfi, as defined in Section [2?T| 
with eigenvalue 1. Due to the odd parity condition T27r,o = —T and, hence, we find that the 
ground state and the proposed state are eigenvectors corresponding to distinct eigenvalues. 

More concretely, it is easy to check that 

(4.1) (V'a,T(27r),^o) = (T2^,oV'a,T(27r),T^o) + ((/-T2^,o)^a,T(27r),^o), 
from which the estimate 

(4.2) K^,,T(2vr),^o)| < ^||(T2.,o-/)^a,T(2vr)|| 
immediately follows. The remainder of this section will be used to prove a bound on 

(4.3) \\Te,oA,T{0)-M^)\\ 
uniformly for 6 G [0, 2tt]. This is the content of Theorem 14.21 

4.1. Observations concerning the twisted ground state. We begin with a warm-up 
exercise involving the twisted ground state. In Section 12.11 we saw that the twisted ground 
state is invariant with respect to the twisted translations; i.e., Tg^_0ipo{6, —6) = 4'o{9, —6*), 
and therefore 

(4.4) de[Tg^-gM0,-0) - M0,-0)] = 0. 

One may rewrite this derivative in the form of an operator acting on '^■^■■> (14.41) is 

equivalent to 

(4.5) D(e)Vo(^,-^) = 0, 
where D[9) is given by 

(4.6) D{e) = deTe,.eTl_s + Te,-e B{e)Tl_, - 3(6). 

Here we have used the differential equation for ipo{6, —0), i.e. (12.201) . and the notation from 
the beginning of Section 13.11 which will be used throughout this section. 

It will be easy to see that the operator D{6) can be written as the sum of two terms, 
Di{9) and D2{6), corresponding to the twists at m and m + L/2, respectively. The goal of 
this subsection is to estimate \\Di{9)ipo{9, —0)\\, see Lemma Wl] below. 

Using (I2.13p . one finds that 

(4.7) (9eTe_eTg _g = i ^ - i ^ Sf^^i^/2+i,v)- 
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One has that D{9) = Di{9) - D2{e) where 

(4.8) D,{e) = I J2 Sfm+i,.) + Te,^eB,ie)Tl_, - 3,(6), 
and 

(4.9) D^id) = I J2 + Te,-eB2{e)Tl_, - B^id), 

For what follows, we will denote by (^A)q = {iPq{6, —9), AiPq{9, —9)) the twisted ground state 
expectation of a local observable A. We have demonstrated in Lemma [2.11 that 

(4.10) = d,Eo{9,-9) = {me,-e)e = {A{9))e, 
for i = 1,2. From this, we conclude that 

(4.11) {To,^oB,i9)Tl^,)^ = {Bm)e = {Md))e = 0, 
as well. Moreover, we similarly have that 

(4.12) (A(^))e = as (5^4,.)) =(E^K-)/ =0' 

for any x G [1,-Z^]. For the last equality above, we used that iPq{9,—9) = W{9)iIjo, W{9) 
commutes with the third component of the spins, rotation invariance implies that the total 
spin is zero, and translation invariance in the 1-direction. 

Since D{9)iJo{9, -9) = 0, we have that Di{9)M0, = D2{9)M9, -9) from which it is 
clear that 

(4.13) = {D{9yD{9))e = 2{Di{9yD,{9))e - 2{Dr{9yD2{9))B- 

As indicated above, we wish to estimate the first term on the right hand side above. We do 
so by estmating the second term. Observe that 

(4.14) {D,{9rDme = - J2 (^(m+i,.) ^(m+L/2+iy))o 

POO 

+ ^ E / (i^M - Sfm+l,v)) {MO))), dt 

+ / / (a,i(Ai(^^))*(a,,(A2(^))-a., (T,%A2(TO,-e))>, c^sdt 

JO 

POO 
poo POO 

+ / / {auiA^{9y*{a^siA2i9y-a,s{Te,-eA2i9)Tl_e)))^dsdt. 
Jo Jo 

That each of these terms is bounded follows from our decay of correlations results found in 
Section [5l3l In fact, we have proven the following lemma. 

Lemma 4.1. Assume Fl, F2, and LSMl-6. If there exists a constant c > such that 
IlL > c for sufficiently large L, then there exists constants C > and k > so that 

(4.15) \\Di{9)M0,-0)f < CL^^^-^e"^^^^, 
for L large enough. 
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Proof. Clearly, one has that 

(4.16) \\D,{e)Mo,-o)\\'= {D,{eyD,{e)), = {D.ieyD^ie)),, 

from (14. 13p above. Applying Theorem 15.81 Lemma [5.91 and Remark 15. 10[ as appropriate, to 
the terms found in (14.141) . one arrives at (14.151) . □ 

4.2. Orthogonality of the excited state. We are now ready to provide the orthogonality 
estimate. 

Theorem 4.2. Assume Fl, F2, and LSMl-6. If there exists a constant c > such that 
IlL > c for sufficiently large L and we choose the parameters a = "^l/ L and T = L/2, then 
there exists constants C > and k > so that 

(4.17) KV^,,T(27r),^o)| < CL^V^^^^ 
when L is large enough. 

Proof. As is clear from (14.21) . Theorem 14.21 follows from bounding the quantity appearing in 
(14. 3 p uniformly for 9 G [0,27r]. A short calculation, using (12.250 . shows that 

(4.18) de[Te,oiJaAO)-MG)] = CaA^) [Tefi^aA^) - ^aA^)] + DaAO)MO), 
where 

(4.19) CaAO) = deTe,o T*^o + ^e,o B^AO) T^o, 
and 

(4.20) DaA0) = deTg^oTg*^, + r,,o 5,,t(^) T,*o, -B^Ad), 

are both anti-Hermitian operators. The first term on the right hand side of (I4.18P is norm- 
preserving, and therefore, we need only bound the norm of the second by Theorem 15. 15[ 

The norm of Da^T{6)'ipa,T{d) will now be estimated by rewriting it in terms of quantities 
for which we have already proven bounds. Each term will be shown to satisfy a bound of 
the form KT7\\ . 

We begin by writing 

\\DaAO)M0)V = TT[DaA0yDaA0)paA0)] 

= TT[DaA0TDaAG)po{o,-e)] 

(4.21) + Tr [DaAOTDaA^) (paA^) - PoiO, -0))] ■ 

The first term on the right hand side above, which is equal to ||Da^j'(^)^/'o(^) ~^)IP; may 
be estimated by comparing it with the vector Di{9)iPq{9, —9) introduced in the previous 
subsection. In fact, 

(4.22) \\DaAO)Md,-d)\\ < pi(^^)^o(^^,-^^)|| + II (DaAd) - D^i9))M^,-0)\\. 

We bounded the first term above in Lemma 14.11 For the second, observe that 

(4.23) DaAG) - Di{9) = Te,-e{BaA0) - Bi{9))Tl^, 

+ Te,o [b,A(^) - B^^^\9))tI, 

-Te-e (^BaA^) " B^^{9)^ Tg_g 
+ B^{9) - BaAO), 
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from which it is clear that 

(4.24) +2||(5,,tW - B^ie)) PH \\. 

That each of these terms satisfies the desired bound follows from Lemma | 
For the final term on the right hand side of fl4.2ip . we insert and remove 
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(4.25) 



iW). 



AW) I 



a local observable supported in A^^\m). Observe that 
(4.26) 



\D 



(W) 
a,T 



where we have used Proposition 15.41 We may write 



(4.27) +Tr \(Da,r{erDaAO) " iPaA^) " Po{0, -0)) 



AW), 



The first term above may be estimated by 



(4.28) 



Tr 



< 



AW) 



AW) I 



where for the final inequality above we used Theorem 13.11 again. 
For the second term, we rewrite the difference as 



(4.29) 

and apply the norm estimate 



(4.30) 

We find that 



(4.31) Tr [D^AOyD^AO) - D^^^^T D':A (0) ) {paA^) - Po{0,-9)) 



< 2 



AW), 



BaAO)-BZ\0) 



< 



4 



BaAd) - B^:^r\e) 



\D„ 



\D 



iW) 



which satisfies the required bound by Lemma 
completes the proof of Theorem 14.21 □ 



and an estimate analogous to (14.261) . This 



5. Auxiliary results 

In this section, we collect a number of auxiliary results, technical estimates as well as a 
few lemmas of a more general nature, which are needed for the proofs in Sections [3] and HI 

We first recall the Lieb-Robinson bounds which are used to demonstrate quasi-locality 
of the dynamics associated to general quantum spin systems, see Theorem 15. 1[ Then, we 
observe in Proposition 15.21 that these Lieb-Robinson bounds may be used to compare the 
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dynamics of a Hamiltonian defined on a given system with the dynamics of the same Hamil- 
tonian restricted to a subsystem. Next, we provide in Lemma [5.31 an exphcit bound which 
apphes to the specific type of interactions we consider in this work. 

In Section [521 we introduce the operators H) which play a prominent role in our 

argument. We first discuss a few of their basic properties, and then use Proposition 15.21 
to estimate the difference that arises in defining the operator with the full Hamiltonian as 
opposed to the Hamiltonian restricted to a subsystems; this is the content of Lemma 15.61 
Lastly, we remark on exactly how this estimate will be used in the main text. 

We review the Exponential Clustering Theorem in Section 15. 3[ and use it to prove a tech- 
nical estimate, see Lemma 15. 9[ Moreover, in this section we also prove Lemma 15.121 This 
result provides an estimate on the quantity || [Ba^ri^, H) — B{A, H)) Po\\ in terms of the 
parameters a, T, and the spectral gap of H, see (15.641) . Here Pq denotes the spectral pro- 
jection onto the ground state of H, and the bound is valid for local observables A satisfying 
PoAPo = 0. 

Lastly, we formulate a statement concerning solutions of certain simple differential equa- 
tions in Section 15.41 

5.1. Lieb-Robinson bounds. For what follows, we adopt the same general framework for 
quantum spin models that was described in Section [L2| including Conditions Fl, F2, and 
the assumption that ||$||a < +oo for some A > (see (II. 7p for the definition of the norm 



We will use the following version of the Lieb-Robinson bound [TB|, which is a variant of 
the results proven in [201 [S]- 

Theorem 5.1 (Lieb-Robinson Bound). Let A > and take $ G B\{A). For any pair of local 
observables A G Ax and B G Ay with X, F C A, one may estimate 

(5.1) \\[c^t{AlB]\\ < ^'l^'jlf" 9x{t) 5^5^F,(d(a:,i/)), 

' x&Xy&Y 



for any t G M. Here {dt] is the dynamics generated by $, and one may take 

[enn^cx{F)\t\ _ ifd{x,Y)>o, 

g2 ii'i'iu Ca{-F) |t| otherwise. 



(5.2) gxit) 



Our proof of the Lieb-Schultz-Mattis theorem relies heavily on comparing the time evo- 
lution corresponding to a given Hamiltonian to that of the Hamiltonian restricted to a 
subsystem. The errors that result from such a comparison can be estimated in terms of a 
specific commutator to which the Lieb-Robinson bounds readily apply. 

We begin with some notation. Let A > and consider $ G B{A). For finite Aq C A, the 
Hamiltonian corresponding to $ restricted to Aq is given by the self-adjoint operator 

(5.3) Ho= J2 "^W- 

XcAo 

We will denote by af'^ the time evolution corresponding to Hq, i.e., for any local observable 
A, aJ°^(A) = e^*^Me-'*-^o for all t G M. 

Proposition 5.2. Let A > and $ G Bx{A.). Suppose the Hamiltonian corresponding to $ 
restricted to a finite volume Aq G A is written as the sum of two self-adjoint operators, i.e., 
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Hq = Hi + H2. Denoting by "* the time evolution corresponding to Hi, for i = 0,1, 2, then 
for any local observable A and t G M, one has that 

(5.4) \\af\A) - ai'\A)\\ < 



/' \\[H2,a'-^\A)]\\ ds. 
Jo 



Proof. Define the function / : M — > ^ by 

(5.5) f{t) := af\A) - a?{A). 

A simple calculation shows that / satisfies the following differential equation: 

(5.6) f\t) =i[h^- Hi, al'^A)] + I [Ho, f{t)] , 

subject to the boundary condition /(O) = 0. As this is a first order equation, the solution 
can be found explicitly: 

(5.7) fit) = (^1^ {t [H„a^^\A)] ) d.^ . 

Using expression (15. 7p and the automorphism property of af'\ it is clear that 

(5.8) l|/(t)|| < /'*' \\[H2,a^^\A)]\\ ds, 

Jo 

as claimed. □ 

To estimate the norm of the commutator appearing in Proposition 15.21 specifically in the 
bound (15. 4p . it is useful to specialize the general Lieb-Robinson bounds described above to 
the exact context we encounter in the present work. For example, we will be interested 
in specific finite volume Hamiltonians, those defined in Section [2] as Hg^gi, and particular 
observables, such as Ai{9) = diHg^^g and A2{0) = d2Hg_g. Let at be the time evolution 

corresponding to the Hg^gi Hamiltonian, and let a[^^ denote the dynamics associated with 
the Hamiltonian H^^} which is defined in (12.81) . We use the following estimate several times. 

Lemma 5.3. Let $ G Bx{A), then there exists constants C > and k > for which 



(5.9) max sup 



*=1'2 6»G[0,27r] 



< Ce'^l*! 



Here it is important that C and k depend only on the properties of the underlying set A and 
the interaction $; they do not depend on the length scale L. 

Proof. We will estimate the above commutator in the case that the observable is Ai{9); an 
analogous result holds for A2{6). Recall that in (12. 9p we wrote as a sum of interaction 
terms. Similarly, if one denotes by Pm{9] Y) := V^{ey(^{Y)Vrn{0) - <^{Y), then Ai{9) may 
be written as 

(5.10) Ai{9) = dePm{0;Y) = J2 Yl ^-(^)* ' 

where K|_ is the set of sites y G F strictly to the right of m. Inserting both of these expressions 
into the right hand side of (15.91) and applying the triangle inequality, it is clear that we must 



(5.11) 
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bound many terms of the form 

$(x), [vuey [si $(r)] v^{9)) 

Term by term, we apply the Lieb- Robinson bound provided by Theorem 15.11 and use that 
the distance between the supports of X and Y is hnear in L; concretely for any x G X and 
y d{x, y) > d{X, Y) > j — 3R. We find that each term described by (15.111) satisfies an 
upper bound of the form 

(5.12) Cit)\MX)\\\Y\\\[Sl,HY)]\\e-''/\ 
where C{t) may be taken as 

(5.13) C{t) 



2\\F\ 



,2Cx{F)\mx\t\+3XR 



Cx{F) 



We need only count the number of terms. The combinatorics of the sums may be naively 
estimated as follows: H^^^ corresponds to a sum of the form 



(5,14) 



< 



E E E + E E E 



n=^-R+l ^^Vl X3{m+n,v) 

whereas for Ai{6) we have that the sum 
(5.15) 



n=^-R+l f eVi X3{m+n,v) 



m+R 



E Es E E E E- 

n=m-R vGVl Y^{n,v) y&Y 



P„(6l;y)^0 



Putting everything together, we have obtained that 



(5.16) 



H 



(S) AW) 



which proves the claim. Recall, 
(5.17) 



< 2C{t) |||$|||i |||$|||2|VL|^(2i?+ l)(2i? 
sup5^||<I>(X)|| 



1$ 



X3x 



and 
(5.1^ 

□ 



$|||2:=supJ]|X| [Sl„<!>{X)] 



x€A 



X^x 



x'€X 



5.2. Approximation of the imaginary time evolution. For our proof of the Lieb- 
Schultz-Mattis Theorem, we introduce an operator which, under certain assumptions, ap- 
proximates the imaginary time evolution corresponding to a given Hamiltonian. In this 
section, we provide several basic estimates of this operator to which we will often refer in 
the main text. 

Let A > 0, $ G B\{A), Aq C A be a finite set, and H be the Hamiltonian corresponding to 
$ restricted to Aq. Denote by at, for t G M, the time evolution determined by H. For any 
local observable A, a > 0, M > 0, and t 7^ 0, define 



(5.19) 



Aa,M{it,H) 



2m 



as{A) -ds, 

M s — It 
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and set Aa{it, H) = limM^oo ^a,M(^^5 H). We use the operator 

(5.20) BaAAH) = - [ Aa{it,H)-Aa{it,Hydt, 

Jo 

to define our variational state in the main text, see (I2.25p . We begin with some basic 
properties. 

Proposition 5.4 (Shanti's Bound). Let $ G B\{A), A be a local observable, a > 0, and 
T > 0. The operator Ba^riA, H) is anti-hermitian and bounded. In particular, 

(5.21) \\BMAH)\\ < ^yf- 

It is important to note that the bound above is independent of the finte volume Aq on 
which the Hamiltonian H is defined. 

Proof. That Ba^riA, H) is anti-hermitian follows immediately from fl5.20p . Combining (15.191) 
and (I5.20p . one finds that 

(5.22) B,,t{A,H) = - r r e-'^^^'^''^aM)^^dsdt, 

^ Jo J-oo S +t 

from which (I5.2ip easily follows as 



oo pT 

>2 



1 , , WAW 



(5.23) \\B^,{AH)\\<^ 1^1 / -^^^^^< , 

□ 

In situations where the local observable A and the Hamiltonian H are fixed, we will often 
write Aa{it) and Bar to simplify notation. The following estimate is a simple consequence 

of dsn. 

Proposition 5.5. Let $ G B\{A) and A be a local observable. One has that 



(5.24) 

Proof. For any t ^ 0, 



T 



Aaiit) - Aa^Miit) dt 



~ 2M x/vra 



(5.25) Aa{it)-Aa^M{it) = %— / a,{A)^—ds, 

2m .l\.\~^M s - it 



2 



and therefore, one has the pointwise estimate 



(5.26) II Aa{tt) - A,^,St) II < e-'^*' e""^^^' 

Upon integration, (15.241) readily follows. □ 

We will now prove an analogue of Proposition 15.21 for the quantities Ba^xi^A, H) introduced 
in (15.201) . The estimate provided below is made explicit in terms of an a priori input, an 
assumed form of the Lieb-Robinson bound, see (15.271) below. 
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Lemma 5.6. Let A > and $ G Bx{A.). Suppose the Hamiltonian corresponding to $ 
restricted to a finite volume Aq C A is written as the sum of two self-adjoint operators, i.e., 
Hq = Hi + H2. Denote by af^ the time evolution corresponding to Hi, for 2 = 0, 1, 2. If, for 
a given local observable A, there exists numbers Ci > 0, i = 1,2,3, for which 



(5.27) 



for all if: G M, then the following estimate holds 



(5.28) 



< cie' 



2T 



C2\t\ - C3 



B,,t{A, Ho) - B,,t{A, Hi) \\ < ^ e""*^ 



A\\ ^ ciM2 



TT 



where M has be chosen as the positive solution of 



(5.29) 



aM^ + C2M - C3 = 0. 



We note that in our applications the numbers q will depend on the observables A and H2; 
in fact, they will be functions of the length scale L. We articulate this dependence explicitly 
in Remark 15.71 below. 



Proof. One may write 

(5.30) Ba,TiA Ho) - BaAA Hi) 



+ 



Aa{lt,Ho)- Aa{lt,Hi)dt 



[ Aa{lt,Hoy - Aa{lt,Hiydt, 

Jo 



and therefore 



(5.31) 



Ba,T{AHo)-Ba,T{A,Hi)\\ < 2 



/ Aa{it,Ho) -Aa{it,Hi)dt 
Jo 



Moreover, the integrand may be expressed as 

Aa{it, Ho) - Aa{it, Hi) = Aa{it, Ho) - Aa,M{it, Ho) + 
(5.32) Aa^Miit, Ho) - Aa,M{it, Hi) + Aa,M{it, Hi) - A^iit, Hi), 

and for j = 0, 1, the bounds 



(5.33) 



[ Aa{lt, Hj) - Aa,M{it, Hj) dt 

Jo 



< 



T \\A\ 



2M 

follow immediately from Proposition 15.51 From this we conclude that for any M > 0, 



\\Ba,T{AHo)-B,,T{A,Hi)\\ < 2 

(5.34) 

Clearly, the pointwise estimate 



/ Aa,M{it, Ho) - Aa,M{it, Hi) dt 
Jo 



2T 



M Jna 



(5.35) 



\Aa,M{it,Ho) - Aa,M{it,Hi)\\ < 



271 



A/ 11^(0) 



ar{A) - ai'\A)\\ 



ds. 



M 
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follows directly from fl5.19p . By Proposition 15.21 we have that 
(5.36) < /"''' \\[H,,a^^\A)]\\ dx, 



JO 

and by assumption fl5.27p . the integrand satisfies a uniform bound for |s| < M. The impli- 
cation is that for all |s| < M, 

\s\ 

Putting everything together, we obtain that 

\\ B^AA Ho) - B^AA H^) \\ < ^^e-^^' 

M \/7ra 



(5.38) + — / / e-'''-'"-''"Usdt. 

71" io J-M 

As M here was arbitrary, we choose it as the (positive) solution of the following quadratic 
equation aM^ + C2M — C3 = 0. In this case, 

(5.39) \\Ba,T[A,HQ) - Ba,T[AHi)\\ < —e 



as claimed. □ 

Remark 5.7. In the main text of the paper, we will use Lemma \5. 6\ for the Hamiltonians 
Ho = Hg _g and Hi = Hg^l each of which depends on the length scale L; see Section 
for the relevant definitions. It is assumed that Hq has a gap 7l > above the ground state 
energy. The local observable A will be exactly as in Lemma \5. 3\ and therefore, the numbers 
Ci, i = 1,2, 3, may be taken as follows: C\ = CL'^^^~^\ 02 = k, and C3 = XL/ A, where again C 
and k depend only on the interaction and the underlying set A. In this case, we will choose 
the parametrization a = 'Jl/L and T = L/2. With this choice, the estimate ( 15. 28\) takes the 
form: 

(5.40) sup B^^rie) - B^^^\d) < CL^'^e"^^-^ (l + ■ 

Here we have used the notation from Section\^ and the fact that the gap 7^ has a uniform 
bound from above; see ( fi.ij) in SectionUi 

5.3. Estimates for gapped systems. We derive two useful results in this subsection. For 
the first we recall the Exponential Clustering Theorem [20119]. and use it to prove a technical 
estimate Lemma [5.91 

The second crucial estimate in this subsection. Lemma 15.121 below, applies specifically 
to gapped systems. It provides a bound on the norm of the difference in the operators 
Ba^T^A, H) and B{A, H) when restricted to the space of ground states corresponding to H. 
The bound applies to local observables A which project off the ground state, i.e. satisfy 
PqAPq = where Pq is the spectral projection of H onto the ground states, and is explicit 
in the parameters a, T, and the spectral gap of H, see (15.641) below. 

We will consider Hamiltonians H, of the type introduced in section ISTTl with an additional 
feautre: a gap above the ground state energy. To state the gap condition precisely, we 
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consider a representation of the system on a Hilbert space Ti. This means that there is a 
representation n : A-^ 13{T-C), and a self-adjoint operator H onH such that 

7r{at{A)) = e**^7r(A)e-**^, 

for all t G M and A & A. For the results which follow, we will assume that H > and that 
i7 G 7i is a normalized ground state, i.e., a vector state for which HQ = and = 1. We 
say that the system has a spectral gap in this representation if there exists 6 > such that 
<j{H) n (0, S) = 0, where cr{H) is the spectrum of the operator H. In that case, the spectral 
gap, 7, is defined to be 

(5.41) 7 = sup{(5 > I a(iJ) n (0,(5) = 0}. 

Let Pq denote the orthogonal projection onto kerif. From now on, we will work in this 
representation and simply write A instead of tt{A). 

The following result concerning exponential clustering was proven in 



Theorem 5.8 (Exponential Clustering). Fix A > 0. Let $ G Bx{A) be an interaction 
for which the corresponding self-adjoint Hamiltonian has a representation H > with a 
normalized ground state vector Q, i.e., HQ = and \\Q\\ = 1. Let Pe denote the family 
of spectral projections corresponding to H . If H has a spectral gap of size 7 > above the 
ground state energy, then there exist fi > such that for any local observables A and B with 
A G Ax, B G Ay, d := dist{X, Y) > 0, and PqBQ = PoB*Q = 0, the estimate 

-fid ( 1 



(5.42) I ( Q, AauiB)Q ) | < C{A, B) e 
holds for allt : < t(4||$|| aCa + 7) < 2\d. Here, one may choose 

(5.43) C{A,B) = \\A\\\\B\\ [ 1 + ^ ^ ^ F (rf(x, y)) + ^ 



and 

(5.44) /i 



4||$||aCa + 7 



7 



The above result easily leads to estimates on integrals of these ground state expectations. 
We state two such bounds in the next lemma, as they will arise in the proof of our main 
result. 

Lemma 5.9. Under the assumptions of Theorem \5.S\ we have the estimates 

(5.45) / \{Q,Aait{B)Q)\ dt < {2fidC{A,B) + \\A\\ ||5||e"^'^) 

Jo 

and 

(5.46) 

-fid 

I {Q,Aai^s+t){B)Q) \dsdt < [{fidf C{A,B) + ||v4|| (2/id + e"'^'^)] — ^. 

Proof. Define T by the equation jT = 2 fid. We have that 

(5.47) / I {Q,Aait{B)Q) \ dt < C{A, B)T e-^"^, 

Jo 
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and also 

(5.48) r I {n,Aau{B)n) \ dt < Mli^g-^^. 

Jt 7 

Combining these two bounds, we arrive at (I5.45p . Similarly, one may estimate 

(5.49) / / I {Q,Aa,^,+t){B)Q) \ dsdt < ^ ' ' e'^^ 

Jo Jo 4 

/.QO /'T/2 WAW llfillT' 

(5.50) / / I {Q,Aa,^,+t)iB)Q) \dsdt < " " " e'^^ 
7t/2 7o ^7 

and finally, 

pco pco WAW II R II 

(5.51) / / I {Q,Aa,^s+t){B)n) \dsdt < " " ' " e"^^. 



T/2 JT/2 7 



□ 



Remark 5.10. In our applications, the Hamiltonian H = Hq _q depends on a length scale 
L and has a gap 7l > above the ground state energy. The support of the observables A 
and B will have a minimal distance d = L/2 — 2R — 1, and moreover, B = B{6) will either 
be Ai{9) or A2{9). In this case, {B{9))g = by Lemma l2J\ and therefore, the assumptions 
of Theorem \5.8{ hold. Here we have used {■)g to denote the ground state expectation corre- 
sponding to ipo{9, —9) It is easy to see that there exist positive constants C and C" for which 
C'IlL < fid < C"L, and thus ultimately constants C and k for which the bounds appearing 
in ( [5.y^5| ) and Ili5.46\ ) niay be estimated by 

(5.52) sup / \{Aait{B{9)))e\dt < C\\A\\\X\—e-^^^^, 

E[0,27r] Jo iL 



9e 
and 

poo poo U^^^ 

(5.53) sup / / \{Aan^s+t){B{9)))e\dsdt < C\\A\\\X\^^ 

6le[0,27r] Jo Jo 7l 



respectively. 

For the next lemma we will need the following basic estimate involving the decay of certain 
Fourier transforms. 

Proposition 5.11. Let a > and T > be given. Define a function Fa^T : —>■ C by 

(6.64) 

For all G M, Fa^riE) > and the estimate 

(5.55) Fa,TiE) < |e-^, 

is valid for E > 0. In the parameter range, E > 2aT > 0, one may also show that 

(5.56) f e-^'dt - FaA-E) < ^e"^. 
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Proof. One may easily verify that for any t > 0, 

(5.57) — / '^ds=^ e-'-e-'^dw, 

2vr2 y_oo s-it 2Vvra 



for all E eM., see e.g. Lemma 1 in [20]. This implies the first claim. Evaluating the Gaussian 
integral yields 

5.58 ^= / e-''"e-^-^dw < -e"", 



in the case that E > 0, from which ( 15.55^ is clear. 

To obtain f l5.56p . we first recall that the Fourier transform of a Gaussian is a Gaussian, 
i.e., for all 2; G C, 



(5.59) e"^ = ^= / e-"^e-"^dx 



and therefore, by rescaling z 1-^ —\^2az, multiplying through by e*^^ (for E G M), and 
changing variables w = ^/2ax + E, we have that 



(5.60) e'^' e-"'' = ^= I e'"^' e-""^ dw 

^ ^ 2v^ ' 



00 



for all z eC. 

Now, by direct substitution into ( ]5.57p . we have that 



(5.61) F,,Ti-E) = —= / e--'" / e-'^e'^dw. 

^^'na Jo 







Applying (15.601) . with the special choice of 2; = it, one sees that 

pT -1 pT /*0 2 

(5.62) / e-^*t;t - F,,t(-^) = / e-"^*' / e"*"" e"^^ du;. 



Jo 2y7ra Jo 

Since w < and t > 0, the integrand above 

(5.63) e e 4a = e e 2" e 4a 

satisfies a trivial bound when E > 2aT. For these values of E, fl5.56p holds. □ 

We may now prove the main estimate for gapped systems. Recall the definitions of the 
operators B = B{A,H) and Ba,T = Ba,T{A,H) in (1X1^]) and ([221, respectively. 

Lemma 5.12. Let H > be a self-adjoint operator and Pe denote the family of spectral 
projections corresponding to H . Suppose H has a gap 7 > 0, and let A be a local observable 
for which PqAPq = 0. If 2aT < 7, then one has that 

(5.64) II (5,^ - B) Po II < Te-t (\\^EA±}^\ + flZ! pp„||. 



7 

Proof. One may rewrite the difference in these operators as 

{Ba,T -B) Po = I {au {A) - Aailt)) dt Po 
Jo 

poo pT 

(5.65) + / aa{A)dtPo+ / Aa{itydtPo, 

Jt Jo 
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Each of these terms may be bounded in norm. 
For any vectors / and g, one may calculate 



(5.66) 



POO noo 

/, / a,t{A)dtP,g) = / {f,e-'" AP,g)dt 

JT JT 



OO fOO 



T J-1 



e-'''d{f,PEAPog)dt, 



where we have used the spectral theorem to rewrite the time evolution and the fact that 
PqAPq = 0. Clearly then, 



(5.67) 

and therefore, 
(5.68) 



/, / au{A)dtPog) < 



ait{A)dtPo 



\APog\ 



e-^' dt 



-7T 



< 



lAPnll. 



7 



Likewise, one may similarly calculate 



T 



(/, / AaittydtPog) 



2m 



{ asiA)f, Pog ) — — ds dt 



s + it 



(5.69) 



Fa,T{E)d{f,PEA*Pog), 



where we have introduced Fa^T{E) = Fa^xiE) as in f l5.54p of Proposition 15.111 above. The 
estimate 



(5.70) 



Aaiitydtp^ 



^ 2 

< -e"^ \\A*Pq\ 
-2 II "I 



readily follows from (15.551) of Propostion 15.111 and the fact that < 7 < -E. 
Lastly, an analogous calculation shows that 



(/, [ait{A) - Aa{it)] Pog)dt 



(5.71) 



-Et 



dt -F, 



-E) 



L "'0 



d{f,PEAPog). 



Thus, for 2aT < 7, we may apply (I5.56P of Proposition 15.111 and establish the bound 



(5.72) 



T 



[au{A) - Aa{it)] dtPo 



< — \\APq\ 
-2 II "I 



Compiling our estimates, we have proven that: if 2aT < 7, then 



(5.73) II {Ba,T - B) Po II < Te-^ 

as claimed. □ 



lAPn 



\A*Pn 



-7T 



7 



\APn 
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Remark 5.13. Applying Lemma l5.12\ to the operator H = Hq ^q, whose spectral projections 
we denote by P^, and the local observable A = Ai{6), we find that there exists a constant 
C > for which, along the parametrization a = 'Jl/L and T = L/2, 

(5.74) sup II iBa,Tie) - B,{e)) PH W <CL^e 

6»e[0,27r] 

5.4. Norm preserving flows. In this section, we collect some basic facts about the solu- 
tions of first order, inhomogeneous differential equations. 

Definition 5.14. Let B be a Banach space. For each ^ G M, let A{6) : B ^ B be a bounded 
linear operator, and denote by X{6) the solution of the differential equation 

(5.75) deX{e) = A{e)X{9) 

with boundary condition X{0) = Xq E B. We say that the family of operators A{9) is 
norm-preserving if the corresponding flow is isometric, i.e., for every Xq G B, the mapping 
■ye : B ^ B which associates Xq X(6'), i.e., 7e(Xo) = X{6), satisfies 

(5.76) II 7e(Xo) II = II Xo II for all 6 eR. 

Two typical examples are the case where i3 is a Hilbert space and A{6) is anti-hermitian 
and the case where ;B is a Banach space of linear operators on a Hilbert space with a spectral 
norm (such as a p— norm with p G [1, +C)o]), and where ^4(6') is a symmetric derivation (e.g., 
i times the commutator with a self-adjoint operator). 

Theorem 5.15. Let A{6), for 6 E R, be a family of norm preserving opeartors in some 
Banach space B. For any bounded measurable function B : R ^ B, the solution of 

(5.77) deY{9) = A{9)Y{9) + B{9), 
with boundary condition Y{0) = Yq, satisfies the bound 

(5.78) II Yi9) - 7,(Fo) || < f \\ B{9') \\ dff . 

Jo 

Proof. For any 6' G M, let X(^) be the solution of 

(5.79) deX{9) = A{9)X{9) 

with boundary condition X(0) = Xq, and let ye be the linear mapping which takes Xq to 
X(^). By variation of constants, the solution of the inhomogeneous equation (15.771) may be 
expressed as 

(5.80) Yi9) = ye (^Yo + ^ {ysyHS (s)) ds^ . 

The estimate fl5.78p follows from fl5.80p as A{9) is norm preserving. □ 

5.5. Existence of local unitaries with vanishing expectation. Consider a finite system 
with a Hamiltonian of the form 

(5.81) F = ^$(X) 

X 

where $ is an interaction as defined at the beginning of the paragraph containing equation 
(11. 7p . In the introduction, (11. ip . we stated a simple upper bound for the spectral gap of any 
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such Hamiltonian with a unique ground state. The argument we gave there made use of a 
one-site unitary U G A{x} with the property that {Q, UQ) = 0. In the following lemma we 
show that such a unitary always exists. 

Lemma 5.16. Let Ti he a complex Hilbert space of dimension at least 2. Then, for any 
density matrix p on Ti, there exists a unitary U G B{l-L) such that TrpU = 0. 

Proof. First consider the case where dimTi is finite and even, or infinite. Let {cq, Ci, . . .} de- 
note an orthonormal basis of eigenvectors of p, with eigenvalues pi ordered in non-increasing 
order. If TC is not separable, it is sufficient that {cq, ei, . . .} contain a basis for the separable 
subspace ranp.) Then, a suitable unitary U can be defined as follows: 

(5.82) [/ = I 

j>0 

If dimH is odd (and hence by our assumptions > 3), it is sufficient to change the first 
summand in (15.821) as follows 

U = a|eo)(eo| - a|ei)(ei| + 6|ei)(eo| + &|eo)(ei| + e"^|e2)(e2| 

+ ® I^2i-i)(e2i| + |e2i)(e2i-i| , 

i>2 

where a, 6 G C and G M, are chosen such that |ap + |6p = 1 and 

e*'^P2 = apo — api . 

This is always possible since p| < (po + Pi)^- 

It is straightforward to check that U thus defined has the desired properties. □ 
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